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OPTIMIZATION  AND  STABILITY  CONTROL  OF  RELATIVISTIC  AND 
PONDEROMOTIVE  SELF-CHANNELING  OF  ULTRA-POWERFUL  LASER 
PULSES  IN  UNDERDENSE  PLASMA 


Introduction 

The  stability  analysis  of  relativistic  and  ponderomotive  self-channeling  of  ultra¬ 
powerful  laser  pulses  in  initially  uniform  underdense  plasmas  has  revealed  [1-4]  three 
major  modes  of  the  self-channeling:  stable  mode,  with  the  formation  of  a  single  stable 
ultra-powerful  channel,  highly  unstable  mode,  which  is  characterized  by  the  development 
of  multiple  peripheral  filaments  with  a  power  in  the  range  of  the  critical  power,  and 
bifurcation  mode  of  the  self-channeling,  when  the  transverse  instability  leads  to  the 
splitting  of  the  beam  into  two  or  more  powerful  channels.  In  the  case  when  the  power  of 
the  laser  beam  is  much  greater  than  the  critical  power  of  the  self-channeling  and  the 
radius  of  the  beam  is  much  bigger  than  the  radius  of  the  channel,  the  highly  unstable 
filamentation  is  the  dominant  mode  of  the  self-channeling.  It  has  been  established 
numerically  and  later  verified  experimentally,  that  unstable  modes  of  the  self-channeling 
can  be  converted  into  the  stable  mode  using  an  appropriate  gradient  of  the  electron 
density  at  the  initial  stage  of  the  channel  formation.  Our  latest  numerical  experiments 
demonstrated  that  this  technique  can  be  used  to  establish  optimization  and  stability 
control  of  multi-TW  and  multi-PW  relativistic  laser  channels  in  underdense  plasma.  The 
stable  multi-PW  laser  channels  contain  the  power  exceeding  104  critical  powers  and  the 
peak  relativistic  channel  intensity  is  in  the  1023  W/cm2  range. 

The  stability  analysis  of  relativistic  and  ponderomotive  self-channeling  is  based 
on  a  relatively  simple  model  [1,5]  that  involves  two  dimensionless  parameters  (r|,  p0) 
representing  the  normalized  power  q  and  normalized  radius  po  of  the  beam  given  by 
rl  =  po/Pcr  and  p0  =  r0cop0/c.  (1) 
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In  equation  (1)  Po  and  ro  denote  the  peak  power  and  the  radius  of  the  incident  beam,  Pcr 
represents  the  critical  power  for  relativistic  and  ponderomotive  self-channeling  given  by 
[1,5] 

©O 

Pc  =  (me20c5  /e2)  jg2(p)pdp  (co/ cop0)2  =1.6198xl010  (to /a>p0)2  W ,  (2) 

0 

where  nie,o,  c,  and  e  have  their  standard  identifications,  go(p)  is  the  Townes  mode  [6],  ©  is 
the  angular  frequency  of  the  laser  radiation,  and  ©p,o  represents  the  angular  frequency  of 
the  laser  radiation,  given  by 

(°P,o  =(4rce2Ne(0/me0)1/2.  (3) 

The  process  of  the  self-channeling  can  be  described  as  the  stabilization  of  the  transverse 
beam  profile  near  one  of  the  z-independent  modes  of  propagation  identified  as  the  lowest 
eigenmodes  [1,5]  Us,o(p)  of  the  governing  nonlinear  Schrodinger  equation.  The  index  s, 
0<s<l,  is  associated  with  normalized  power  r|  of  the  eigenmodes,  p  =  rcoPio/c  represents 
the  normalized  transverse  variable,  where  r  is  the  transverse  coordinate.  The  normalized 
radius  pe,o  of  the  eigenmodes  Us,o(p)  can  be  defined  as 


P,.»»[2]uJo(p)pdp/U!1„(0)]1'2.  (4) 

0 

Since  eigenmodes  Us,o(p)  represent  the  normalized  transverse  field  profile  of  the 
relativistic  channels,  the  eigenmode  curve  pe,o(n)  identifies  the  normalized  radius  of  the 
relativistic  channel  as  a  function  of  normalized  power  trapped  in  the  channel.  Thus,  the 
radius  of  the  channel  rCh  can  be  expressed  as 


and 


CPe,0 

®r 


(5) 


2 


(6) 


Since  pe,0  varies  slowly  [1]  for  r\  >  L5,  we  conclude  that  the  radius  of  the  relativistic 
channel  scales  as 


rch  ~  Ne 


-1/2 


(7) 


From  the  equations  (4)  and  (6)  the  peak  laser  intensity  in  the  relativistic  channel  Ich  can 
be  expressed  as 

4k  N„ 


^"p2  NePch^2’ 

Ke,0  iycT 


(8) 


where  PCh  represents  the  power  trapped  in  the  channel.  Again,  since  pe,o  varies  slowly  [1] 
for  r)  >  1.5,  the  peak  intensity  in  relativistic  channels  scales  as 


Ich  ~  ^  Pch  XT2. 
N„ 


.(9) 


For  a  constant  ratio  of  Ne/Ncr  we  have  from  equations  (6)  and  (9)  the  simple  scaling 
relations 

rCh  ~  X, ,  Ich  ~  Pch  ^  •  (10) 

The  relativistic  and  ponderomotive  self-channeling  of  ultra-powerful  (t|  =  P0  /  Pcr  »  1) 
laser  pulses  with  initial  radius  of  the  beam  r0  much  bigger  than  the  radius  of  the  channel 
(ro  »  rCh,  or  po  »  pe,o),  is  typically  highly  unstable  and  results  in  the  formation  of 
multiple  peripheral  filaments.  In  order  to  establish  stable  channel  formation  and  optimize 
the  efficiency  of  the  power  compression,  the  good  match  between  the  incident  radius  of 
the  beam  and  the  radius  of  the  formed  channel  has  to  be  arranged: 


ro  ~  rCh  or  po  ~  pe,o- 


(11) 
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The  condition  (11)  can  be  achieved  dynamically  by  the  adjustment  of  the  local  electron 
density  Ne  (since  rCh  ~  Ne'1/2)  with  an  appropriate  longitudinal  gradient  in  the  electron 
density. 

The  major  parameters  of  the  conducted  numerical  stability  analysis  are:  the 
incident  peak  power  and  the  radius  of  the  laser  beam,  the  level  of  azimuthal  perturbations 
in  the  initial  transverse  intensity  distribution,  the  length  of  the  initial  stage  of  the  self- 
channeling  with  longitudinal  gradient  in  the  electron  density,  the  level  of  electron  density 
variation  and  the  shape  of  electron  density  profile  at  the  initial  stage,  where  stable  ultra¬ 
powerful  channels  are  formed.  The  stability  analysis  was  performed  for  the  perturbed 
Gaussian  beams  with  incident  transverse  amplitude  distribution  given  by 

U0(r, cJj)  =  I"2  exp(-0.5(r/r0)2  X(l  +  (r /r0)4£e9  cos (#))•  (12) 

The  level  of  weak  azimuthal  perturbations  in  equation  (12)  is  characterized  by  the 
parameter  5pert  defined  as 

dpert  =  max  1 10  (r,  ft  )  - 10  (r,  <p2  )  |  / 10 .  (13) 

For  sq  =  0.01  (q  =  1-4)  in  equation  (12)  we  have  5pert  =  0.063,  and  for  eq  =  0.05  (q  =  1-4) 
the  level  of  azimuthal  perturbations  is  8pert  =  0.393. 

Results 

Figure  1  illustrates  the  highly  unstable  mode  of  relativistic  and  ponderomotive 
self-channeling  in  initially  uniform  plasma  with  the  formation  of  multiple  peripheral 
filaments.  The  initial  transverse  intensity  distribution  is  presented  in  Figure  1(a).  The 
level  of  azimuthal  perturbations  in  incident  transverse  intensity  profile  is  5pert  =  0.063. 
The  major  laser  beam  and  plasma  parameters  in  this  case  are:  X  =  248  nm,  Po  =  39  TW,  ro 
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=  10  Jim,  Xo  =  1.2  x  1019  W/cm2,  sq  =  0.01  (q  =  1-4),  Spert  =  0.063,  Ne,0  =  3.8  x  102°  cm'3. 
The  normalized  power  and  beam  radius  in  this  case  are  r\  =  P0  /  Pcr  =  50,  po  =  20.  Figure 
1(b)  depicts  the  transverse  profile  with  multiple  peripheral  relativistic  filaments  evolved 
in  the  process  of  highly  unstable  relativistic  and  ponderomotive  self-channeling.  The 
stability  of  the  self-channeling  and  efficient  power  compression  can  be  established  by 
introduction  of  multi-stage  self-chanfieling.  The  first  stage,  which  initiates  the  stable 
channel  formation,  involves  the  exponential  longitudinal  gradient  in  electron  density 
profile  defined  by  Ne,0(z)  =  Ne>0(0)exp(otz)  with  Ne,o(0)  =  1.9  x  1019  cm’3,  Ne,o,max  =  3.8  x 
10  cm  .  All  parameters  of  the  laser  beam  are  the  same  as  in  the  previous  case  presented 
in  Figure  1.  The  dynamics  of  stable  self  channeling  of  the  laser  beam  from  Figure  1(a)  in 
the  plasma  with  longitudinal  electron  density  profile  is  presented  in  Figure  2.  It  is  evident 
that  highly  unstable  structure  with  multiple  filaments  (see  Figure  1(b))  has  been 
transformed  into  the  formation  of  stable  ultra-powerful  laser  channel.  Our  numerical 
modeling  demonstrated  that  similar  stable  channel  formation  can  be  established  in  the 
multi-PW  range.  Figure  3  illustrates  the  formation  of  stable  multi-PW  relativistic  laser 
channel  at  the  initial  stage  of  relativistic  and  ponderomotive  self-channeling  involving 
longitudinal  gradient  in  the  electron  density.  The  major  laser  beam  and  plasma 
parameters  in  this  case  are:  X  =  248  nm,  P0  =  7.8  PW,  r0  =  2  pm,  Io  =  6.2  x  1022  W/cm2,  eq 
=  0.01  (q  =  1-4),  8pert  =  0.063,  Ne>0(0)  =  5.0  x  1019  cm'3,  Ne,0,max  =  3.8  x  102°  cm’3.  The 
initial  power  of  the  laser  pulse  Po  is  Po  =  104  Pcr,  and  the  power  trapped  in  the  stable 
relativistic  channel  PCh  is  in  the  same  range.  In  this  case,  the  longitudinal  profile  of  the 
electron  density  at  the  initial  stage  of  the  relativistic  and  ponderomotive  self-channeling 
has  an  exponential  shape  defined  by  Ne,0(z)  =  Ne,0(0)  x  exp(ocz).  The  performed 
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numerical  experiments  demonstrated,  that  for  laser  pulse  and  plasma  parameters 
corresponding  to  Figure  3,  stability  of  the  multi-PW  relativistic  channel  formation  can  be 
established  with  the  various  shapes  of  the  electron  density  profile  at  the  initial  stage  of 
the  self-channeling.  One  of  the  results  is  presented  in  Figure  4,  which  demonstrates  the 
process  of  stable  relativistic  and  ponderomotive  self-channeling  with  electron  density 
profile  defined  by 

Ne,o(z)  =  Ne,o(0){  Ne,o,max/Ne,o(0)  —  (Ne,o,max/Ne,o(0)  — l)Xexp[-(z/Zexp,2)  ]}>  (14) 

with  Ne>o,max/Ne,o(0)  =  7.6  and  zexp,2  -  300  pm.  All  major  laser  pulse  and  plasma 
parameters  are  the  same  as  in  the  case  presented  in  Figure  3.  Stable  multi-PW  relativistic 
laser  channel  presented  in  Figure  4  contains  power  in  the  order  of  104  critical  powers. 

The  process  of  stable  formation  of  multi-PW  relativistic  channels  with  the  power 
exceeding  104  Pcr  is  presented  in  Figure  5.  The  major  laser  beam  and  plasma  parameters 
in  this  case  are:  X  =  248  am,  P0  =  9  PW,  r0  =  5  pm,  Io  =  1.1  x  1022  W/cm2,  eq  =  0.01  (q  = 
1-4),  5Pert  =  0.063,  Ne,o(0)  =  7.6  x  1018  cm'3,  Ne,o,max  =  3.8  x  1020  cm'3.  The  longitudinal 
profile  of  the  electron  density  at  the  initial  stage  of  the  relativistic  and  ponderomotive 
self-channeling  has  an  exponential  shape  defined  by  Ne,o(z)  =  Ne,o(0)  x  exp(az).  The 
transverse  profile  of  the  incident  laser  beam  and  transverse  profile  of  evolved  stable 
multi-PW  relativistic  laser  channel  are  presented  in  Figure  5(a)  and  Figure  5(b), 
respectively.  Figure  5(c)  illustrates  the  process  of  stable  channel  formation.  The 
trajectory  of  stable  relativistic  and  ponderomotive  self-channeling  in  the  plane  of 
dimensionless  parameters  (t|,  po)  (see  equation  (1))  is  presented  in  Figure  5(d).  It  is 
evident  from  the  position  of  point  Bstabie  that  the  stable  multi-PW  relativistic  channel 
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depicted  in  Figure  5(b)  contains  more  than  104  critical  powers.  Note,  that  the  peak 
intensity  in  the  channel  exceeds  1023  W/cm2. 

Our  numerical  experiments  demonstrated  that  the  same  approach,  involving  initial 
stage  of  the  self-channeling  with  longitudinal  electron  density  profile,  results  in  stable 
multi-PW  relativistic  channel  formation  even  in  the  case  of  severely  distorted  incident 
laser  beams  with  relatively  big  initial  beam  radius.  Figure  6  illustrates  one  of  these 
examples.  The  major  laser  beam  and  plasma  parameters  in  this  case  are:  X  =  248  nm,  P0  = 
9  PW,  r0  =  5  pm,  Io  =  1.0  x  1022  W/cm2,  eq  =  0.05  (q  =  1-4),  5pert  =  0.393,  Ne,o(0)  =  9.5  x 
1017  cm 3,  Ne>o,max  =  3.8  x  IO20  cm"3.  The  longitudinal  profile  of  the  electron  density  at  the 
initial  stage  of  the  relativistic  and  ponderomotive  self-channeling  has  an  exponential 
shape  defined  by  Ne,o(z)  =  Ne>o(0)  x  exp(az).  The  strongly  azimuthally  perturbed 
transverse  intensity  profile  of  the  incident  laser  beam  is  presented  in  Figure  6(a).  The 
transverse  profile  of  the  evolved  stable  multi-PW  relativistic  channel  is  depicted  in  Figure 
6(b).  Figure  6(c)  demonstrates  the  dynamics  of  stable  channel  formation  in  the  case  of 
severely  azimuthally  aberrated  incident  multi-PW  laser  beams  with  relatively  big  initial 
beam  radius. 

Summary 

In  conclusion,  detailed  studies  of  the  transverse  stability  of  relativistic  and 
ponderomotive  self-channeling  have  demonstrated  that  multi-stage  self-channeling,  with 
the  initial  stage  involving  an  appropriate  gradient  in  the  electron  density  profile,  results  in 
stable  formation  of  multi-PW  relativistic  channels  in  underdense  plasma  with  power 
exceeding  104  critical  powers  and  peak  channel  intensity  in  the  1023  W/cm2  range. 
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CAPTIONS 


Figure  1.  Highly  unstable  relativistic  and  ponderomotive  self-channeling  of  azimuthally 
perturbed  multi-TW  laser  pulse  in  initially  uniform  underdense  plasma.  The  laser  and 
plasma  parameters  are:  X  =  248  nm,  P0  =  39  TW,  r0  =  10  pm,  Io  =  1.2  x  1019  W/cm2,  Ne>0 
=  3.8  x  IO20  cm'3.  The  normalized  power  and  beam  radius  in  this  case  are  T|  =  Po  /  Pcr  = 
50,  po  =  20.  The  incident  azimuthally  perturbed  transverse  amplitude  distribution  is 
defined  by  the  equation  (12)  with  eq  =  0.01  (q  =  1-4).  The  corresponding  level  of  weak 
azimuthal  perturbations  is  characterized  by  8Pert  from  the  equation  (13),  8per,  =  0.063. 
Panel  (a)  shows  weekly  azimuthally  perturbed  transverse  intensity  profile  of  the  incident 
laser  beam.  Panel  (b)  demonstrates  the  transverse  structure  with  multiple  peripheral 
filaments  evolved  in  the  process  of  highly  unstable  relativistic  and  ponderomotive  self¬ 
channeling. 

Figure  2.  Dynamics  of  stable  multi-TW  relativistic  channel  formation  at  the  initial  stage 
of  the  self-channeling  with  exponential  longitudinal  electron  density  profile  defined  by 
Ne>o(z)  =  Ne,0(0)exp(az)  with  Ne,o(0)  =  1.9  x  1019  cm'3,  Ne,0>max  =  3.8  x  IO20  cm'3.  All  laser 
pulse  parameters  are  the  same  as  in  the  previous  case  depicted  in  Figure  1:  X  =  248  nm, 
Po  =  39  TW,  ro  =  10  pm,  Io  =  1.2  X  1019  W/cm2.  The  incident  azimuthally  perturbed 
transverse  intensity  profile  is  the  same  as  the  distribution  from  Figure  1(a)  with  sq  =  0.01 
(q  =  1-4)  and  8pert  =  0.063  (see  equations  (12)  and  (13)).  The  dynamics  of  stable 
relativistic  channel  formation  is  presented  in  (x,z)  plane,  where  z  is  the  axis  of 
propagation  and  x  is  one  of  the  transverse  coordinates. 
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Figure  3.  Dynamics  of  stable  multi-PW  relativistic  channel  formation  at  the  initial  stage 
of  the  self-channeling  with  exponential  longitudinal  electron  density  profile  defined  by 
Ne,o(z)  =  Ne>o(0)exp(az)  with  Ne,o(0)  =  5  x  1019  cm'3,  Ne>o,max  =  3.8  x  1020  cm'3.  Laser 
pulse  parameters  are:  X  =  248  nm,  P0  =  7.8  PW,  r0  =  2  pm,  Io  =  6.2  x  1022  W/cm2.  The 
incident  azimuthally  perturbed  transverse  amplitude  distribution  is  defined  by  the 
equation  (12)  with  eq  =  0.01  (q  -  1-4).  The  corresponding  level  of  azimuthal 
perturbations  is  5pert  =  0.063  (see  equation  (13)).  The  power  trapped  in  the  stable 
relativistic  multi-PW  laser  channel  is  in  the  range  of  104  critical  powers. 

Figure  4.  Dynamics  of  stable  multi-PW  relativistic  channel  formation  at  the  initial  stage 
of  the  self-channeling  with  longitudinal  electron  density  profile  defined  by  the  equation 
(14)  with  Ne>o(0)  =  5  x  1019  cm'3,  Ne,o,max  =  3.8  x  102°  cm'3,  and  zexp,2  =  300  pm.  Laser 
pulse  parameters  are  the  same  as  in  the  previous  case  illustrated  in  Figure  3:  X  =  248  nm, 
Po  =  7.8  PW,  ro  =  2  pm,  Io  =  6.2  x  1022  W/cm2.  The  incident  azimuthally  perturbed 
transverse  amplitude  distribution  is  defined  by  the  equation  (12)  with  sq  =  0.01  (q  =  1-4). 
The  corresponding  level  of  azimuthal  perturbations  is  5pert  =  0.063  (see  equation  (13)). 
The  power  trapped  in  the  stable  relativistic  multi-PW  laser  channel  is  in  the  range  of  104 
critical  powers. 

Figure  5.  Stable  multi-PW  relativistic  channel  formation  at  the  initial  stage  of  the  self¬ 
channeling  with  exponential  longitudinal  electron  density  profile  defined  by  Ne,o(z)  = 
Ne,o(0)exp(az)  with  Ne,o(0)  =  1.9  x  1019  cm'3,  Ne>o,max  =  3.8  x  102°  cm'3.  Laser  pulse 
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parameters  are:  X  =  248  nm,  P0  =  9  PW,  r0  =  5  pm,  I0  =  1.1  x  1022  W/cm2.  The  incident 
azimuthally  perturbed  transverse  amplitude  distribution  is  defined  by  the  equation  (12) 
with  sq  =  0.01  (q  =  1-4).  The  corresponding  level  of  azimuthal  perturbations  is  5pert  = 
0.063  (see  equation  (13)).  Panel  (a)  shows  the  incident  azimuthally  perturbed  transverse 
intensity  profile.  The  transverse  profile  of  the  evolved  stable  multi-PW  relativistic 
channel  is  presented  in  Panel  (b).  Panel  (c)  demonstrates  the  dynamics  of  stable  channel 
formation  in  (x,z)  plane,  where  z  is  the  axis  of  propagation  and  x  is  one  of  the  transverse 
coordinates.  The  trajectory  of  stable  relativistic  and  ponderomotive  self-channeling  in  the 
Cn>  Po)  plane  (see  equation  (1))  is  depicted  in  Panel  (d).  The  power  of  the  evolved  stable 
multi-PW  relativistic  channel  exceeds  104  critical  powers.  The  peak  channel  intensity  is 
in  the  1023  W/cm2  range. 

Figure  6.  Stable  relativistic  channel  formation  in  the  case  of  severely  azimuthally 
aberrated  incident  multi-PW  laser  beam  with  relatively  big  initial  beam  radius.  The 
exponential  longitudinal  electron  density  profile  is  defined  by  Ne,o(z)  =  Ne,o(0)exp(az) 
with  Ne>o(0)  =  9.5  x  1017  cm'3,  Ne,o,max  —  3.8  x  1020  cm"3.  Laser  pulse  parameters  are:  X  = 
248  nm,  Po  =  9  PW,  ro  =  5  pm,  Io  =  1.0  x  1022  W/cm2.  The  incident  severely  distorted 
transverse  amplitude  distribution  is  defined  by  the  equation  (12)  with  eq  =  0.05  (q  =  1-4). 
The  corresponding  level  of  azimuthal  perturbations  is  8pert  =  0.393  (see  equation  (13)). 
The  strongly  azimuthally  perturbed  transverse  intensity  profile  of  the  incident  laser  beam 
is  presented  in  Panel  (a).  The  transverse  profile  of  the  evolved  stable  multi-PW 
relativistic  channel  is  presented  in  Panel  (b).  Panel  (c)  demonstrates  in  (x,z)  plane  the 
dynamics  of  stable  relativistic  and  ponderomotive  self-channeling  of  severely 
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azimuthally  aberrated  incident  multi-PW  laser  beam  with  relatively  big  initial  beam 
radius.  The  power  of  the  evolved  stable  multi-PW  relativistic  channel  exceeds  104  critical 
powers.  The  peak  channel  intensity  is  in  the  1023  W/cm2  range. 
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Figure  1(a) 


8  nm,  P()  =  39  TW,  r0  -  1 0.0  jum,  P 
2  x  1019  W/cm2,  Spcrtlirbalions  =  0.063 
3.8  x  102()  cnr3 
ly  Uniform  Electron  Density 


P0  =  7.8  PW 
r0  =  2.0  |im 
l0  =  6.2  x  1 022  W/cm2 
Ne0(O)  =  5  x  I019  cnr3 
Ne;o.,„?x  =  3.8xl02»cm-3 
Longitudinal  Gradient  in 
Electron  Density 


Channel 


SitaiMe  j.ceiaEOTStne  met  cttod'eromoiMve 
Ctemelleg  of  Asmotoaiy  Ferterlbed  MoMATW 
Laser  Fuilses  m  IJedlerdeese  Plasma 


_  Ai. 


fTN  }} 


A.  =  248  nm 

P0  =  7.8  PW 

r0  =  2.0  |i.m 

I0  =  6.2  x  1022  W/cm2 

Ne0(0)  =  5  x  1019  cur3 

NeAn,ax  =  3.Sxl02«Cm-2 

Longitudinal  Gradient  in 
Electron  Density 


Channel 


Nc,o(z  )  -  Nc,o,max  [l-(l-NCio(0)/NcAmax)*exp(-(z/zcxpJ2)  ] 

NcFOyNcAmax  =  1/7-6,  Zc.\p,2  =  300  hm 


Figure  4 


17 


blmbm  IJMafewnste  ma.4.  .^oiudleromoflve  Self- 
Gfemdimg  oir  AuBjioJitUhi^Ily  I^'erterbed  Multi-TVS 
Laser  Pnlseo  m  JEiderdleiase  Plasma 


leaser  if 


A,  =  248  nm,  P0 
I0=  1.1  x  1022 

Ne„(0)  =  7.6  x 


9.0  PW,  r0  =  5.« 

/rrn  2  S 

7  ’  u perturbations 

ils  cur3,  Ne  0  niax 


Longitudinal  Gradient  in  Electro 
Ne.o(z)  =  Ne  0(0)  x  exp(az) 
Ne,o,maX  Ne-0(0)  =  50,  zmax=  150( 


ntICE! 


*/; 


Figure  5(b) 
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Figure  6(c) 


